Time-dependent fields of a current-carrying wire 
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Abstract. The electric and magnetic fields of an infinite straight wire carrying 
a steady current which is turned on abruptly are determined using Jefimenko's 
equations, starting from the standard assumption that the wire is electrically 
neutral in its rest frame. Some nontrivial aspects of the solution are discussed in 
detail. 



1. Introduction 

Consider an infinite straight linear wire carrying the current I(t) = for t < 0, 
and I(t) = Iq for t > 0. That is, a constant current Io is turned on abruptly at 
time t — 0. What are the resulting electric and magnetic fields? This apparently 
simple electrodynamic problem is posed and solved as Example I 0.2 in the excellent 
textbook of Griffiths pQ. Starting from the standard assumption that the wire is 
electrically neutral in its rest frame, without or with the current, which implies that 
the scalar potential V is zero, the retarded vector potential A is calculated, and 
then the electric and magnetic fields are obtained according to E = —dA/dt and 
B = V x A, respectively. (Note that in this case the Coulomb and Lorenz gauges 
lead to same potentials since the wire is electrically neutral.) While this solution 
is correct, we believe that the problem has some intriguing aspects and as such 
deserves further attention. In this note we present a solution to the problem using 
Jefimenko's equations and point out some pitfalls which could be dangerous for novices. 
Thus, hopefully, our analysis should be instructive for advanced undergraduate and 
beginning graduate students. 

2. Solution using retarded potentials 

For the convenience of the reader, we first give the solution using retarded potentials, 
in some more detail than that given in Griffiths's book. 

As is well known, the retarded vector potential A(r,t) at field point r and time 
t is given by 

f J(r',t r ) 3 , 



A(r,*) = p / , x ' dV, (I) 
Air J \r - r'\ 

where J{r',t r ) is the current density at a source point r' and the retarded time 
t r = t — \r — r'\/c. Let the infinitely long wire lie along the z axis. The current in the 
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wire, which is assumed to have an infinitesimal cross section, is turned on abruptly at 
t = 0, and thus the current density can be expressed as 

j(M) = f^e(t)i, (2) 

2n s 

where s is the distance from the wire, 6(s) is the one-dimensional Dirac delta function 
normalized as L 5(s) ds = 1 and 0{t) is the Heaviside step function, 

«*>-{!: 2!^ < 3 > 

The setup is not realistic, but, in principle, it could be realized approximately with 
a large superconducting loop of negligible cross section in an inhomogeneous axially 
symmetric magnetic field, the symmetry axis coinciding with the axis of the loop. 
If the loop, initially at rest with no current, is moved quickly along the symmetry 
axis into a new resting position, a persistent current is produced in it, since the total 
magnetic flux through the loop is constant (cf, e.g., [2]). 

Equations (1) and (2) imply that the vector potential A at a distance s from the 
wire is given by 



47T Jn ,/n ./_™ 27T s' d 



voio . [°° eit-V^+z^/c) 



z 



. , dz'. (4) 

4tt ./-oo Vs 2 + z' 2 

Here, cylindrical coordinates s, tp, z are used and d=[s 2 +s' 2 — 2ss'cos((j>— <fi')+(z— z') 2 ]^ , 
which is the distance between the field point (s, 0, z) and a source point (s', 0', z'); in 
the second line, the integration with respect to s' and a transformation z — z' — > z' 
reduce the distance to (s 2 + z' 2 ) 1 / 2 . As demanded by the problem's symmetry, the 
vector potential is independent of z and <j>. The Heaviside-function factor in the 
integrand of the integral in the second line of (4) causes the potential to vanish at 
times t < s/c and limits the integration interval to the values of z' satisfying 



\z'\ < yjc 2 t 2 - s 2 . (5) 
Integral (4) for the vector potential A thus evaluates aa3 



A( s ,t) = ^zQ(t- s /c) \ 



o 



[\n(ct + y/c 2 t 2 - s 2 ) - \ns]@(t - s/c). (6) 



2vr 

The electric field is therefore given by 
E{s,t) = -dA(s,t)/dt 

= _Moio c e(t _ s/c) (?) 

2vr ^c 2 t 2 - s 2 y ' ' W 

X The standard convention is understood according to which f(x)@(x — xq) = whenever x < a;o, 
even when the expression /(x) happens not to be defined at these values of x. 
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and the magnetic field bj[| 



B(s,t) = V x A(s,t) = -(dA z /ds)d> 



2ns 



0- 



ct 



rO(t-s/ C ). 



(8) 



In both (7) and (8), the delta- function terms that arose from the derivatives of the 
Heaviside step function in (6) dropped out on account of the property 

f(x)S(x-xo) = f(x )S(x-xo) (9) 

of the delta function. Inspecting equations (7) and (8) we see that the fields E(s,t) 
and B(s,t) attain in the limit t — > oo their familiar static values and ([IqIq/2'ks)4>, 
respectively, and that both these fields diverge when t — > s/c. 



3. Solution using Jefimenko's equations 

As is now well known, starting from the retarded solution to the inhomogeneous wave 
equations for the fields E and B [H [5] , or from the familiar retarded potentials [TJ [S] , 
or otherwise [7J, the time-dependent generalizations of the Coulomb and Biot-Savart 
laws can be derived: 



E(r,t) 



B(r,t) 



1 



47re n 



Mo 
4tt 



g(r',U) 



e(r',t r ) n J(r',t r ) 



cTZ 2 



c 2 TZ 



J(r',t r ) J(r',t r 



TV 



cTZ 2 



xKd 3 r', 



(10) 



(11) 



where g is the volume charge density, TZ = r — r' , and the dots denote partial 
differentiation with respect to time. These equations, showing explicitly true sources 
of E and B, were first derived by Jefimenko |3]. We now shall calculate the fields E 
and B in the problem at hand using Jefimenko's equations. 
In our problem, the charge density vanishes, 

g = 0, (12) 

since by assumption the wire is electrically neutral, and using equation (2) we get 



2ir s 

A 8 ') 



J(s',t r ) = ^e(t-s/cy- 



5(t - Vs 2 + z' 2 /c) z. 



(13) 



(14) 



Here, cylindrical coordinates are used again, and the delta-function property (9) and 
the taking, with no loss of generality, the field coordinate z to be reduced the 
retarded time t r to the same value as that in equation (4); the step- function factor 
in (14) expresses the fact that the step function in (13) entails that not only the 

§ It is perhaps worthwhile to note that expression (6) for the vector potential A(s, t) resembles the 
quasi-static vector potential at a distance s from the midpoint of a straight wire of finite length 
21 = 2\/c 2 t 2 — s 2 carrying a constant current 7o- -A calculation of the magnetic field according to 
B = V X A in which the distance-dependent length 21 is treated as a constant would be equivalent 
to the use of the Biot-Savart law. However, this law is not applicable beyond the quasi-static regime 
(cf., e.g., [3]). 
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current density itself but also its partial time derivative vanishes for times t < s/c. 
Substitution into Jefimenko's equation (10) then gives 

47TC0C 2 Jo S ' J_ 00 Vs 2 + z' 2 

which can be evaluated easily using the decomposition of the delta function [5J [5] 

S(t - Vs 2 + z' 2 /c) = cH \5(z' - \Jc 2 t 2 - s 2 ) + 5{z' + yjc 2 t 2 - s 2 )} (16) 

\/c 2 t 2 — s 2 L J 

as 

E(s, t) = - - C 0(f - s/c), (17) 

in full agreement with the electric field (7), obtained using the retarded vector 
potential. In a similar fashion, using equations (11), (13), (14) and (16) we obtain 

u( a ^ n« dz ' x 1 r 8(t-VF+**/c) 

B(s,t) = — — &(t—s/c)s(p / — , + - / t, d, 

4tt [J-V^t^ (s 2 + z l2 f/ 2 c J_ 00 s 2 + z' 2 

2ns v y/c 2 t 2 - s 2 V ' ' ^ ' 

in full agreement with the magnetic field (8), obtained using the retarded vector 
potential. 

4. Discussion 

At first sight, the fact that the fields E and B obtained diverge when t —> s/c+ while 
vanishing for t < s/c may seem disturbing. However, a closer examination reveals 
that, in the correct solution to the problem, E and B must tend to infinity when 
t — > s/c+. This is more transparent through the use of Jefimenko's equations (10) 
and (11) than the use of a retarded vector potential. It is clear from equations (13), 
(14) and (16) that the abrupt turning of the current on at t = necessarily yields 
an infinite time derivative of the current density, producing in the fields a cylindrical 
'shock wave' that diverges at the time t = s/c at a distance s from the wire. Similar 
to the instructive example of Jackson of an abruptly turned on electric dipole [9] , the 
diverging fields E and B are here artifacts of the unphysical, instantaneous turn-on 
of the current (jj 

The divergences disappear if the current is not turned on abruptly but is increased 
gradually during a short time interval t. As a simple example, assume that the current 
increases linearly from zero at t — to a steady non-zero value Iq at t > r, replacing 
accordingly expression (2) for the current density by 



J M = o 

2ir s 



T T 



(19) 



|| We remind the reader that a similar situation is found in the well-known iiC-circuit problem 
of charging a capacitor of capacitance C by connecting it instantaneously to a constant voltage V 
through a resistor of resistance R, assuming that the charge Q on the positive plate is zero at t = 0. 
The standard (tacit) assumption that the inductance L of the circuit is zero then leads to the equation 
V = QIC + RI. In the unphysical setup of the problem (the abrupt closing of a circuit with L = 0), 
the correct solution must satisfy the unphysical initial condition I(t = 0) = V/R, despite the fact 
that the current / vanishes for t < 0. 
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Using Jefimenko's equations, the resulting fields are then obtained to be 
Et(s, t) = - [l n (ct/s + ^c 2 t 2 -s 2 /s) 6(t - s/c) 

Z7T TV. \ J 

- In (c(t -t)/s + Vc 2 (f-r) 2 -s 2 /s) 0(i - t - s/c) (20) 

and 

B T (a, i) = ^£ ^_ [Vc 2 t 2 - s 2 6(t - s/c) - V^ft-r) 2 - s 2 G(t - t - s/c)} . (21) 
2irc st L 



v 2 a(m) - 4,^^ - -^^e(*) *. (22) 



While the fields ([20)1 and (|2Tj) are finite for any non-zero parameter r, their limits 
t — y can be shown easily to be the fields (7) and (8), respectively, that diverge when 
t -> s/c. 

There is another query. The retarded vector potential and Jefimenko's equations 
are both derived under the assumption that the sources (charges and currents) are 
localized in a finite region of space, but our problem involves an infinitely long current- 
carrying wire. Therefore, the question arises as to the validity of the solution foundQj] 
However, inspecting equation (6) we see that for any given finite s and t only a 
finite segment of the wire contributes to the retarded vector potential. Figuratively 
speaking, retardation makes the infinitely long wire finite. Still, one should check that 
the vector potential (6) satisfies the requisite inhomogeneous wave equation, 

1 d 2 A(s,t) = mq/q 8(a) , 
dt 2 2tt 

and that the fields (7) and (8) satisfy all Maxwell's equations. 

Let us check first whether Maxwell's equations are satisfied. Straightforward 
calculations yield that the fields (7) and (8) are divergenceless, confirming the 
equations 'V-E = q/cq, where g—0, and V-.B = 0. While a straightforward calculation 
of the curl of the electric field (7) confirms that the fields obtained satisfy Faraday's 
law, a similar calculation of the curl of the magnetic field (8) using the standard 
cylindrical-coordinate formula 

VxF( S )0=^£ (23) 
s os 

appears to yield only that VxB = eo/iodE/dt, instead of the full Ampere-Maxwell 
law, which reads in our case 

Vx B = ^^e W £ + Wo f. (24) 
Z7T s at 

Here, however, it is important to bear in mind that when the sources of a field are 

idealized point, line or surface distributions of charge and/or current, described by 

generalized functions such as the Dirac delta function, great care must be taken 

to employ in the field or potential differential equations generalized (distributional) 

derivatives instead of the usual (classical) oneo Keeping this in mind, a careful 

examination of the differential operation on the magnetic field (8) implied by formula 

1 Recall that in electrostatics the standard solution to the Poisson equation is not generally valid for 
charge distributions extending to infinity (cf., e.g., |10l 111] ). 

+ For example, calculating the Laplacian of the potential 1/r of a unit point charge using classical 
derivatives yields V 2 (l/r) = for all r > 0; at r = 0, the classical Laplacian is simply not defined. In 
the well-known relation V 2 (l/r) = —Aiz8(r), the Laplacian is in fact the generalized one, as it must be 
since this relation expresses equality of two generalized functions. To avoid confusion, some authors 
denote generalized differential operators by a bar, writing thus V 2 (l/r) = — 47r<5(r) 12. 1131 IT4] . 
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(|23p reveals that it involves an expression that can be written as the Laplacian of the 
natural logarithm of the cylindrical coordinate s0 



ld_ 

s ds 



ld_ 

s ds 



s-Qns) 



V 2 (lns). 



(25) 



In terms of the usual (classical) derivatives, this Laplacian vanishes for all s > and 
is not defined at s = 0, but in view of the fact that the current density involves the 
delta function we should use here the relation 



i\\ns)= S M 



(26) 



Employing this relation in the evaluation of V x B, it is found easily that the fields 
(7) and (8) now satisfy the full Ampere-Maxwell law (f24)). 

Since the relation (|26p may appear to be novel to some readers, we give an informal 
proof of it using a limiting procedure in which In s is regularized as In y/s 2 + a 2 and the 
limit a — > is taken after integrating the product of V 2 In Vs 2 + a 2 and a well-behaved 
'test' function /(s): 



lim 

a->0 



/•oo 

/ V 2 ln Vs 2 + a 2 f{s)sds 
Jo 



lim 

a->0 



/•OO 


' 2a 2 " 


la 


_(s 2 + a 2 ) 2 



f(s)sds. 



(27) 



The expression in square brackets is the Laplacian of In \J s 2 + a 2 , which is now a well- 
behaved function of s for any a ^ 0; its integral over the whole plane is independent 
of a, equalling 2n. Splitting the integral into integrals over intervals < s < S and 
S < s < oo so that the function /(s) can be expanded in the first interval in a Taylor 
series around s = 0, the limit (l27l) can be evaluated as 



lim 



a-»0 



2a 2 



2 



E 



= lim 

a-tO 



+ a 2 ) 2 ^ 

' 71—0 

/(0) 



/ (n) (0)f 

71 



s ds + lim 



a-s-0 



2a 1 



f(s) 



(s 2 + a 2 ) 2 



( ds 



S 

a arct S- - £2 



Sa 2 



f'{0) + O(a 2 ,a 2 lna) + 



= /(0). (28) 



_l+a 2 /S* 2 

The second limit in the first line vanished since the integral of f(s)s/(s 2 + a 2 ) 2 , where 
/(s) is assumed to be integrable over the whole interval (0, oo), is guaranteed to 
converge to a value that remains finite when a — > 0. Having thus shown that 



lim 

o^0 



/•oo 

/ V^nVs 2 +a 2 f{s)sds = f(0), 
Jo 



where /(s) may be any well-behaved function of s, we can write 
lim V 2 In yjs 2 + a 2 = S(s) /s, 



(29) 



(30) 



which is the relation ^ with V 2 (lns) = lim a ^ V 2 In Vs 2 +a 2 (Note that the 



relation (|26l) may be interpreted as the Poisson equation for the electrostatic potential 
— (l/27reo) Ins of a charge density S(s)/2irs, which is that of an infinitely long straight 
line of charge of unit line charge density QUI EJ US].) 

Checking whether the vector potential (6) satisfies the inhomogeneous wave 
equation (I22[) is a somewhat cumbersome but in every step straightforward calculation. 

* This fact is perhaps more transparent when using VxB = V(V • ^4) — V 2 A and recognizing that 
for A expressed by equation (6), V 2 (lns) appears explicitly in the calculation of V 2 A. 
U Strictly speaking, the limit here and in l|30|l is the weak limit (cf., e.g., [8]) and the Laplacian in 
(1 26 I t is the generalized (distributional) one (cf., e.g., |12l ). 
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Using the relation (1261) and making extensive use of relation (9) in the terms containing 
S(t — s/c), which arise through differentiations of 0(t — s/c), we obtain 

V 2 A( S , t) - -L^l = -^ 6 Me(t - s/c) z, (31) 
c A ot z 2ir s 

which confirms, as it musi0, that the vector potential (6) satisfies the wave equation 

(|22j) since, on account of (9), 5(s)Q(t-s/c) = 6(s)G(t). 

Closing the discussion, we remark that the solution to the problem of finding the 

fields due to an infinite straight linear wire carrying a constant current that is turned 

on abruptly appears to be relevant to the related problem of finding charges and fields 

in a current-carrying wire of finite cross section, where it appears that an infinite time 

is needed for establishing a stationary charge distribution [15 . However, it should be 

noted that the time-dependent fields obtained here approach their values at t — > oo 

very rapidly. For example, the field (18) at a distance s = 1 m and a time t = 3 

/zs differs from its asymptotic value by less than 1 part in a million (the wire is then 

required to be at least 2 km long). 

5. Conclusions 

Wc calculated the electric and magnetic fields of an infinitely long wire carrying 
a constant current that is turned on abruptly, using Jefimenko's equations. Our 
calculations confirmed the terse solution given to this problem in an example of 
Griffiths 's text, but our method appears to provide more insight than the standard 
approach via the retarded potentials: the divergence of the resulting fields E(s,t) 
and B(s,t) when t — > s/c then can be seen easily as a necessary consequence of the 
'unphysical' setup of the problem. We also calculated the fields for a more realistic 
case in which the current in the wire is increased linearly in a nonzero time to its 
final steady value. We believe that our analysis was an instructive demonstration of 
not only the power but also the possible pitfalls of using the delta function in the 
calculations of electrodynamics. 

Acknowledgments 

We thank David Griffiths for illuminating correspondence and the anonymous referee 
for constructive comments. DVR acknowledges support of the Ministry of Science and 
Education of the Republic of Serbia, project No. 171028. VH co-authored this paper 
in his private capacity; no official support or endorsement by the Centers for Disease 
Control and Prevention is intended or should be inferred. 

References 

[1] Griffiths D J 1999 Introduction to Electrodynamics 3rd edn (Upper Saddle River, NJ: Prentice 
Hall) 

[2] Romer R H 1990 The motion of a superconducting loop in an homogeneous magnetic field: the 

harmonic oscillator equation in an unfamiliar setting Eur. J. Phys. 11 103-6 
[3] Hnizdo V 2003 Comment on 'About the magnetic field of a finite wire' Eur. J. Phys. 24 L15— L16 
[4] Jefimenko O D 1989 Electricity and Magnetism 2nd edn (Star City, WV: Electret Scientific) 
[5] Jackson J D 1999 Classical Electrodynamics 3rd edn (New York: Wiley) 

ft The facts that V x B = [i J + t fi dE/dt, B = V x A, E = -dA/dt and V • A = already 
imply that CIA = —fj,oJ. 



Time- dependent fields of a current- carrying wire 



8 



[6] Griffiths D J and Hcald M A 1991 Time-dependent generalizations of the Biot-Savart and 

Coulomb laws Am. J. Phys. 59 111-7 
[7] Ton T C 1991 On the time-dependent, generalized Coulomb and Biot-Savart laws Am. J. Phys. 

59 520-8 

[8] Vladimirov V S 1979 Generalized Functions in Mathematical Physics (Moscow: Mir) 

[9] Jackson J D 2010 Comment on 'Maxwell equations and the redundant gauge degree of freedom' 

Eur. J. Phys. 31 L79-L84 
[10] Redzic D V 2012 The calculation of the electrostatic potential of infinite charge distributions 

Eur. J. Phys. 33 941-6 

[11] Jimenez J L, Campos I and Roa-Neri J A E 2012 Confusing aspects in the calculation of the 

electrostatic potential function of an infinite line of charge Eur. J. Phys. 33 467-71 
[12] Kanwal R P 2004 Generalized Functions: Theory and Applications 3rd cdn (Boston: Brikhauser) 
[13] Estrada R and Kanwal R P 1995 The appearance of nonclassical terms in the analysis of point- 
source fields Am. J. Phys. 63 278-78 
[14] Hnizdo V 2011 Generalized second-order partial derivatives of 1/r Eur. J. Phys. 32 287-97 
[15] Redzic D V 2012 Charges and fields in a current-carrying wire Eur. J. Phys. 33 513-23 
[16] Stratton J A 1941 Electromagnetic Theory (New York: McGraw-Hill) 



